Is the Z+(4430) a radially excited state of D s l 
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We present the interpretation that the recently discovered Z + (4430) by the Belle Collaboration 
can be a radial excitation of the cs state, being consistent with an observed value of the product of 
branching ratios, B(B° -> K*Z ± (4430)) x £(2* (4430) -> ir^ij)') ~ 10" 5 . We give an explicit cs 
candidate for this state by calculating the mass value in our semirelativistic quark potential model 
and also give a natural understanding for the facts that the decay mode Z — > J /ifjix + has not yet 
been seen while Z — > ip'n can be seen. 

PACS numbers: 12.39.Hg, 12.39.Pn, 12.40.Yx, 14.40.Lb, 14.40.Nd 
Keywords: potential model; spectroscopy; heavy mesons 

A series of exotic X, Y, and Z charmonium-like mesons have been discovered by the B factories, among which 
the recent discovery of Z + (4430) by the Belle Collaboration [l[ draws attention of many physicists because of the 
following reasons; Z + (4430) is charged and hence cannot be a cc charmonium state or a ccg hybrid meson and it might 
£N) . be the first charged tetraquark state because it is too heavy to be a charged Qq meson and furthermore, strangely 
enough, only the decay mode Z + — > i//7r + is found while the mode Z + — > J/ipir + has not yet been seen. The product 
of branching ratios is determined to be 

B{B° -> K T Z ± {U30)) x B(Z ± (U30) -► ttV) = (4.1 ± l.O(stat) ± 1.4(syst)) x 10" 5 . (1) 

I/"") ■ Mass and width of this particle are measured as 

o : 

mz = (4433 ± 4 ± 2) MeV, T z = ( 45 +™ + J MeV. (2) 



Looking at the above data, many people regard this particle Z + (4430) as a strong and plausible candidate of a 
tetraquark state. According to such a point of view, the second factor of branching ratios given by Eq. ([I]) is in the 
0(1), because the decay width of Z + (4430) is 45 MeV as given by Eq. © and the decay ^±(4430) -> tt^' occurs 
through strong interaction, and then the first factor becomes the 0(1O -5 ) from Eq. {I), whose suppression is due to 
a very small recombination probability for making four quarks (ccqq) into a tetraquark state. There already appear 
several papers on this state whether a tetraquark or molecular state interpretation is possible or not. However, 

all models presented so far cannot give a convincing explanation on why the mode Z + — > J/i/jit + has not been seen, 
though some people [H partly account for the suppression of Z + — > J /vbir + based on the D\D* or D[D* resonance 
model of Z + (4430). Furthermore, it should be noted that the decay width given by Eq. (JTJ is not a partial decay 
width for the process Z + — > J/tpir + but a total decay width of Z + (4430), because it is obtained from fitting the 
invariant mass distribution of ir + + ip' to the Breit-Wigner resonance formula. 

In this letter, we would like to propose a different interpretation from these papers, i.e., this state can be a higher 
radial excitation of D s state. In this model, the first factor, B(B° — > K T Df (4430)) , of the branching ratios of Eq. (JT} 
replacing Z ± (4430) with D^(4430) is of the order 0(1) because the numerator contains a gluon interaction as shown 
in Fig.l below, while the second factor, B(Df (4430) — > ^ip'), is of the order 0(1O -5 ) because an excited £>+(4430) 
can decay to D°K + /D + K° and also to their excited particles D^K^ via strong interactions with the decay width 
of tens of MeV because of its large mass value 4430 MeV, in addition to the decay (4430) — > ip' + ir + occurring via 
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TABLE I: First order mass spectra of the radial excited D s (n = 2, 3, 4) in I/ttlq. Units are in MeV. 
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TABLE II: Optimal values of parameters. 

Parameters a™= 2 ' 3 ' 4 a (GeV^ 1 ) b (GeV) m u , d (GeV) m s (GeV) m c (GeV) m b (GeV) 
0.344 1.939 0.0749 0.0112 0.0929 1.032 4.639 



weak interactions. Based on our semirelativistic quark potential model pdl|-[l6l|. we give not only the numerical mass 
value corresponding to Z + (4430) but also give the reason why the decay mode Z + — > J/i[;ir + has not yet been seen. 

The masses of radially excited D s (cs) states are calculated and presented in Table U by using our semirelativistic 
quark potential model, which succeeds in reproducing the mass levels of all existent heavy mesons including recently 
discovered higher states of D, D s , B and B s [H, Ei]. In Table HI the same physical parameters as in Refs. [la . [Til . fl5j | 
are used and given in Table HT1 except for a™ =1 where i = 1,2, 3, and 4 denote t he p rincipal quantum number. In this 
paper we have adopted the same value a\ for i = 3 and 4 as that for i = 2 [l4[ . From Table HI one may identify 
one of two states, n = 4 3 Pi(l+) 4440 MeV, and n = 4 3 P 2 (2+) 4411 MeV, as a candidate for Z+(4430). Here we 
have discarded the last two columns of Table U because these states, 3 -Di(l~) and " 3 Z)2(2~)" in each row, are largely 
affected by the off-diagonal elements of the interaction Hamiltonian among these states with the quantum number 
2S+1 Lj in the larger state space 11]. Furthermore, it is natural to consider that the Z + (4430) does not have a large 
orbital angular momentum but that it is rather the S state or at most the P state, even if it has large n. 

The parity is not a good quantum number in the decay process D s (4430) — * iJj't: + since this process goes via a weak 
interaction as shown in Fig. 3(a) and furthermore the spin of the excited D s is 1 if the particles ?/>' and n + originated 
from the excited D s are in the relative S state, while it can be 2, 1 or if they are in the relative P state. Thus, 
both of n — 4 " 3 Pi"(l + ) and n — 4 3 P2(2 + ) can be a possible candidate. This should be compared with the case of 
Z + (4430) to be a tetraquark, in which a decay of tetraquark to ip'+ir occurs through a strong interaction and hence 
the allowed spin-parity of a tetraquark should be 1 + for the relative S state of i[>' and ?r + or 2~, 1~, and - for the 
relative P state of tp' and n + because the parity is conserved in a strong interaction process. 

Next, if the Z + (4430) is a radially excited state of D s , how can we explain the suppression of the mode Z + — > J/tpn? 
The answer is as follows. Let us consider the decay process D+(4430) — > ip'ir or Z?+(4430) — > J/ipn as (cs) — * (cc) +ir 
by treating the tt as an elementary Nambu-Goldstone boson. Then the decay amplitude is proportional to the 
overlapping integral of the wave functions of (cs) and (cc) states. We notice that if the node of an initial cs state 
wave function is the same as that of a final cc state wave function, then the decay amplitude is expected to be large. 
On the other hand, if they are different, the magnitude of the decay amplitude is small. Here, we assume that the 
initial state is a radially excited D s state with higher node. Then, if the final state is J r /ip whose node is zero, the 
magnitude of the decay amplitude is small or negligible. To see how it works, let us assume the trial wave functions 
for D s , J/ip, and "0' being expressed by the Hermite polynomials, 'I'x^y, z), as 

* x (x, y, z) = ^-H 2m {m x x)H 2m {m x y)H 2m (m x z)e -™ 2 A* 2 +y 2 +* 2 )/2 ^ (3) 

where H2 m {x) is the 2m-th Hermite polynomial with a node m and we have confined the wave function in the range 
0<^<oo,0<y<oo, and < z < oo. Nq is a normalization and mx is mass of the particle X, which is introduced 
to make the arguments of the Hermite polynomials dimensionless. The ratio of decay rates of D S (A430) — > J/tpir to 
D s (4430) -> ip'ir is proportional to | / d 3 x^ Xm ^ xo\ 2 /\ J <Px^ X m^xi\ 2 where m, 0, and 1 are nodes of the D s (4430), 
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FIG. 1: Decay process through D 3 : B° -> Z + (U30)K 
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FIG. 2: Decay process through tetraquark: B° -> ^+(4430)^" -> ip'^K' 



J/f/;, and "0': respectively. The calculated result is roughly given by 



r(D s (4430) J/^tt) 
r(D s (4430) -> jP'tt) 



6.7 x 1(T 4 
8.5 x 1(T 6 



(n 
(n 



6.9 x 1(T 7 (n = 4) 



= 2 for D s ) 
= 3) 



(4) 



where n is the principal quantum number of the initial _D S (4430) state and its node number is given by ti — 1. We 
know that the node of J/tp is zero and that of ip' is one, hence if we assume the node of D s (4430) is three with n = 4, 
then, the ratio of the production rate for the mode Z — > J/ipn to Z — ► becomes 6.9 x 10~ 7 , which is negligibly 
small. 

Some comments are in order for distinguishing the cases of Z = tetraquark and Z = radially excited D s . The chain 
decay process B° — > Z + (4430)-ft'~ — > ip'ir + K~ for Z + (4430) being a radially excited D s or a tetraquark goes through 
the Feynman diagrams shown in Fig. 1 and Fig. 2, respectively. As mentioned early, the excited Z? s (4430) decays via 
strong interactions as D s (4430) -> DK (D+(4430) -> D + k /D°K + ). In addition to this mode, Z? s (4430) -> D^K^ 
are also possible where Z?^*- 1 and JsT^*) are excited states of D and K, respectively, if they are kinematically allowed. 
These channels will be identified with D + K production accompanying one or more pions and they are also signal 
channels for the Z + (4430). However, among these channels, Z? s (4430) — > DK will be dominant because of the phase 
space effect. 

1. Figure 1 in the case of Z = radially excited D s involves two weak bosons, hence this is the second order process 
in the weak interaction. On the other hand, Fig. 2 in the case of Z = tetraquark involves only one weak boson 
and thus this is the first order weak interaction process. However the latter case needs to be multiplied with a 
recombination probability to form a tetraquark, which must be O(10 -5 ) as mentioned early. Thus, it is expected 
that these two cases would be in the same order of magnitude for a chain decay rate. Therefore, we cannot 
distinguish a tetraquark and a radially excited D s as long as we look at the product of branching ratios alone 
given by Eq. ([T]). 

2. As one can easily notice from Fig. 2, three charged states, ± and 0, of Z(4430) must be available if Z is a 
tetraquark state because [ccud\, [ccdu], [cc(uu — dd)/y/2], and [ce(uu + dd) /V2] are possible. On the other hand, 
if Z is a radial excitation of D s , then only the ± charged states are possible because only cs and cs are possible 
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FIG. 3: D s Decay 




FIG. 4: Tetraquark Decay 



states. Namely if Z is a tetraquark state, its isospin is one and/or singlet, while if Z is a radially excited state 
of D s , then its isospin is 1/2. 

3. To distinguish these two possibilities, i.e., the case of a tetraquark and the case of a radially excited D s , let us 
consider the diagrams, Figs. 3 and 4, for Z = £> s (4430) and Z = T(4430) (T means a tetraquark), respectively. If 
we assume Z = D s (4430), then the process D+(4430) — > ip'n + (Fig. 3(a)) goes through one weak boson while the 
processes £>+(4430) -> D a K+(Fig. 3(c)) / D + K° (Fig. 3(d)) are strong decays. The process D+(4430) -> D a D+ 
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(Fig. 3(b)) needs to go through one weak boson exchange, hence the following ratios are expected to be obtained. 

r(L>+(4430) -> D + K°/D°K + ) T (L>+(4430) -> D°D+) 

r(D+(4430) ->iP'tt+) > ' r(Z?+(4430) -> iP'tt+) ~ ^ ' 

On the other hand, if we assume Z = T(4430), either the process T+(4430) -> D°K+ (Fig. 4(c)) or T+(4430) -> 
D + K a (Fig. 4(d)) goes through one weak boson exchange. The process T + (4430) — ► D°D + (Fig. 4(b)) is not 
a weak interaction but a strong interaction process. Hence the following ratios between these processes and the 
process T + (4430) — > V' 71 " 4 " (Fig. 4(a)) are expected to be obtained. 

r(T+(4430) -> D + K°/D°K + ) T (T+(4430) — > Z)°D + ) 

r(T+(4430) -v 7//7T+) <1 ' r(T+(4430) ^ ~ 1 ' 

Therefore if one measures the process Z — > Z?if and its ratios to Z-> ^V, i.e., Eqs. |(5J) and (jSJ), then one can 
distinguish these two possibilities whether Z is £> s (4430) or T(4430). 

4. Since a radially excited D s decays to tp' + n via weak interactions, the partial decay width r(D+(4430) — ► ip'ir) 
would be very small, contrary to the one for Z+(4430) = T+(4430), r(T+(4430) — > ip'n), which is via a strong 
interaction decay. However, as shown in Eq. ((SJ, Z)+(4430) will dominantly decay to D + K° / D° K + and hence 
the total decay width of (4430) becomes an order of a strong decay. As described in the beginning, in the 
Belle experiment only the total decay width of Eq. @ being in the order of strong interactions, is obtained 
by fitting the invariant mass of tp' + 7r to the S'-wave Breit-Wigner resonance formula. Therefore, the present 
experiment does not rule out the possibility of Z + (4430) being the Z?+(4430). 

The Belle Collaboration discovered the rather narrow resonance Z + (4430) in the invariant mass distribution of 
tp' + 7r + in the decay mode B° — > K~i/j'tt + . However, in the B° decay there are other decay modes, B° — > 
K-K°D + /K-K + D° or B° -> R-D D + , in addition to this special mode B° -► A'-Z+(4430) -> K~^'tt+. If 
the Z + (4430) is a radially excited D s state, one can find Z + (4430) more frequently in the invariant distribution of 
K°D + /K + D° in the mode B° -> K~ K a D+J K~ K+ D a . On the other hand, if the Z+(4430) is a tetraquark, it 
must be very difficult to find it in the mode B° — > K~ K° D + / K~ K + D° . Therefore, this type of decay mode is very 
important to determine whether Z + (4430) is a tetraquark or a radially excited D s . 

Recently the CDF Collaboration reported the observation of the decay mode B^ — * J / ipn^ |Ti} ■ This must become 
a very interesting example to test our model: in this decay mode the heavy mesons in the initial and final states have 
both zero nodes and therefore, this decay rate should be very large compared with the rate for the process B^r — * V ; ' 7r± 
with the node of tp' to be 1, which might be strongly suppressed or not be observed. 

In this letter, we have discussed the possible interpretation of Z+(4430) as a radially excited D s state. We presented 
the mass levels of radially excited D s states and gave the reasonable explanation on why the mode Z — > J/ ipn was 
not observed in the Belle experiment. As shown in Table |H we can see a wealth of excited D s states. We would 
like to stress that the experimental search for resonances in DK and £)WifW invariant mass distributions is very 
interesting not only for observing Z + (4430) but also for discovery of n = 2 and n = 3 excited D s states. Related to 
this, it is remarkable that the n = 2 3 5i(l^) 2715 MeV and n = 2 3 P (0 + ) 2856 MeV were already observed and 
reproduced well by our semirelativistic model, as shown in Table 1. 

We would like to urge the analysis on the invariant mass distribution of K° D + / K + D° in the decay mode B° — > 
K~ K° D + 1 K~ K + D° which must contain a fruitful physics of excited D s states. 
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